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Turbulent stress is the fundamental quantity in the filtered equation for large-scale veloc- 
ity that reflects its interactions with small-scale velocity modes. We develop an expansion 
of the turbulent stress tensor into a double series of contributions from different scales of 
motion and different orders of space-derivatives of velocity, a Multi-Scale Gradient (MSG) 
expansion. We compare our method with a somewhat similar expansion — due to Yeo and 
Bedford, and to Leonard — that is based instead on defiltering. Our MSG expansion is 
proved to converge to the exact stress, as a consequence of the locality of cascade both in 
scale and in space. Simple estimates show, however, that the convergence rate may be slow 
for the expansion in spatial gradients of very small scales. Therefore, we develop an ap- 
proximate expansion, based upon an assumption that similar or 'coherent' contributions 
to turbulent stress are obtained from disjoint subgrid regions. This Coherent-Subregions 
Approximation (CSA) yields an MSG expansion that can be proved to converge rapidly 
at all scales and is hopefully still reasonably accurate. As an important first application 
of our methods, we consider the cascades of energy and helicity in three-dimensional tur- 
■ bulence. To first order in velocity-gradients, the stress has three contributions: a tensile 

stress along principal directions of strain, a contractile stress along vortex lines, and a 
shear stress proportional to 'skew-strain.' While vortex-stretching plays the major role 
in energy cascade, there is a second, less scale-local contribution from 'skew-strain'. For 
helicity cascade the situation is reversed, and it arises scale-locally from 'skew-strain' 
while the stress along vortex-lines gives a secondary, less scale-local contribution. These 
conclusions are illustrated with simple exact solutions of 3D Euler equations. In the first, 
^ , energy cascade occurs by Taylor's mechanism of stretching and spin-up of small-scale 

vortices due to large-scale strain. In the second, helicity cascade occurs by 'twisting' of 
rS \ small-scale vortex filaments due to a large-scale screw. 



1. Introduction 

It is well-recognized that turbulent cascades are essentially multi-scale phenomena, 
and involve the coupling of modes at distinct scales. The property of locality implies 
that these interactions are mainly between adjacent scales (Eyink (2005)). Of course, 
this locality is rather weak and modes at scales differing by an order of magnitude from 
a fixed scale can make a substantial contribution to transfer across that scale. In the 
filtering approach (Germano (1992)) the nonlinear interaction between scales is embodied 
in the turbulent stress tensor r that appears in the equation for the velocity u, low-pass 
filtered at length-scale t. This stress is the contribution to spatial transport of large- 
scale momentum generated by quadratic self-coupling of the subfilter scales. Because of 
the locality property, the subfilter modes that contribute most to the stress are those 
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at scales only somewhat smaller than £. This fact raises the hope that a constitutive 
relation may be constructed for the stress, if the adjacent small scales can be somehow 
estimated from the resolved modes. Unfortunately, it is not hard to see that eliminating 
the small scales produces a stress whose dependence upon the large-scale velocity u is, 
in general, spatially nonlocal, history-dependent, and stochastic (Lindenberg, West & 
Kottalam (1987), Eyink (1996)). Thus, an exact constitutive relation for the turbulent 
stress is formally available but it is quite unwieldy and not of direct practical utility. 

In apparent contradiction to these remarks, several authors have developed a closed 
constitutive formula for the turbulent stress (Yeo & Bedford (1988), Leonard (1997), 
Carati, Winckclmans & Jeanmart (2001)). The basic idea of their approach is to defilter 
u to obtain the unfiltered velocity field u and then to use the latter to calculate the 
stress. For many common filter kernels, it is possible to evaluate the resulting formula 
as a concrete expansion in powers of the filtered velocity-gradients. Although nothing is 
proved about the convergence of this series, it seems in principle to provide a solution to 
the 'closure problem' of turbulence. However, as we argue at length below, this solution is 
illusory for several reasons. Most obviously, defiltering is not defined for all filter kernels. 
Furthermore, the defiltering operator, even when defined, is unbounded on the natural 
function spaces for the velocity field (e.g. finite-energy functions). Thus, the convergence 
cannot hold in general. 

Nevertheless, it is an extremely attractive idea to develop an expansion for the stress in 
powers of the filtered velocity-gradients. Similar expansions have proved useful in many 
areas of physics, e.g. for one-particle distribution functions in the solution of the Boltz- 
mann equation (Enskog (1917), Chapman & Cowling (1939)) or for Ginzburg-Landau 
free energies of superconductors (Gorkov (1959), Tewordt (1965)). We shall here develop 
a convergent gradient expansion for the turbulent stress. It is somewhat more intricate 
than the expansion developed in Yeo & Bedford (1988), Leonard (1997) and Carati et 
al. (2001), since it is expressed by a summation simultaneously over the order of space 
gradients and over an integer index indicating the scale of motion involved. Thus, it is a 
multi-scale gradient expansion. This series expansion will be proved below to converge, 
as a consequence of the locality of the turbulent cascade both in space and in scale. 
Of course, the rate of convergence may be slow, especially for the Taylor expansion in 
space of small-scales, so that very high order gradients could be required to obtain an 
accurate result. We diagnose the reasons for potentially poor convergence, and, on that 
basis, develop also an approximate expansion which will converge rapidly at all scales. 
This approximation may give reasonable accuracy with just a few low-order gradients. 

Because it is multi-scale, the expansion considered here does not by itself give a closed 
'constitutive' relation for the stress. However, it may be a useful point of departure in 
developing a closure for the stress, if supplemented with a scheme to estimate subfilter 
velocity-gradients in terms of filtered velocity-gradients. This in line with some large- 
eddy simulation (LES) approaches which construct subgrid-stress models by creating 
'surrogate' subgrid modes. See Domaradzki & Saiki (1997), Misra & Pullin (1997), Scotti 
& Meneveau (1999), Burton et al. (2002), and, for an extensive review, Domaradzki 
& Adams (2002). Our emphasis in this paper is on fundamental physics rather than 
on closure models, but we hope to pursue this in future work. Even without a closure 
prescription, the formula we develop for the stress makes many testable predictions. 
Concrete conclusions will be deduced here for the joint cascade of energy and helicity 
in three space dimensions (3D) and, in a following paper [Eyink (submitted)], for the 
inverse energy cascade predicted in two dimensions by Kraichnan (1967). 

The contents of this paper are as follows: In Section 2 we briefly review the filtering ap- 
proach to turbulence. In the main Section 3 we develop our Multi-Scale Gradient (MSG) 



Multi-scale Gradient Expansion 



3 



expansion for the turbulent stress. In Section 4 we develop a more rapidly convergent 
but less systematic approximation, which we call the Coherent-Subregions Approximate 
Multi-Scale Gradient (CSA-MSG) expansion. In Section 5 we present the application of 
our method to 3D energy and helicity cascades. Technical proofs and calculations are 
given finally in four Appendices. 



2. Filtering Approach and Turbulent Stress 

We first give a general discussion of the mechanics of energy transfer between scales in 
a turbulent flow. Following Germano (1992), we resolve turbulent fields simultaneously 
in space and in scale using a simple filtering approach. We consider initially an arbitrary 
dimension d of space. Thus, we define a low-pass filtered velocity 

u(x) = J d d r G e (r)u(x + r), (2.1) 

where G is a smooth mollifier or filtering function, nonnegative, spatially well-localized, 
with unit integral J d d r G(r) = 1. The function Gi is rescaled with £, as Gg(r) = 
£~ d G(r/£). Likewise, we can define a complementary high-pass filter by 

u'(x) = u(x)-u(x). (2.2) 

If the above filtering operation is applied to the incompressible Navier-Stokes equation 

d t u+ (u-V)u = -Vp + v A u, (2.3) 

with V-it = determining the pressure p, then one obtains 

a f M+(^V)tI=-V'T-Vp + i/AiI, (2.4) 

where 

T ii/ii uli (2-5) 

is the stress tensor from the scales < I removed by the filtering. 

The equation for energy balance in the large scales is (Piomelli et al. (1991), Eyink 
(1995)): 

<9 t e + V-J = -II-is\Vu\ 2 , (2.6) 

with large-scale energy density e = (1/2)U 2 , spatial energy transport vector in the large 
scales J = (e + p)u + u-t — i/Ve, and scale-to-scale energy flux 

II = -Vu:t. (2.7) 

The latter quantity is the rate of work done by the large-scale velocity gradient against 
the small-scale stress. Of course, it may be rewritten in various equivalent forms as 

n = -Vu:t = -S:r = -S:r. (2.8) 

The first follows from incompressibility of the velocity field, where 

T = T-(Trr)l/d (2.9) 

is the so-called deviatoric stress (with / the d x d identity matrix). The second follows 
from symmetry of the stress tensor, where 

is the large-scale strain rate. The third follows from both properties combined. 
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The following formula 

r = J d d rG e (r)Su(r)Su(r)- J d d r G e (r)6u(r) ■ J d d r G e (r)Su(r). (2.11) 

represents the stress as a tensor product of velocity increments Su(r; x) = u(x+r) — u(x) 
averaged over the separation vector r with respect to filter function Ge(r) at length-scale 
I. It is easily verified by multiplying out the increments and integrating (Constantin,E & 
Titi (1994), Eyink (1995)). This expression implies, as a direct consequence, the matrix 
positivity of the stress (Vreman, Geurts & Kuerten (1994)). It was also the crucial point of 
departure in our discussion of scale locality properties in Eyink (2005). This same formula 
shall play a central role in our development of the multi-scale gradient expansion in this 
work. A decomposition of (2.11) that we shall find useful is 

t = q-u'v! (2.12) 

where 



and 



q{x) = J d d rG e {r)Su(r;x)Su(r;x), (2.13) 

u'(x) = - J d d r G e {r)Su(r; x). (2.14) 

It is easy to check that u'(x) in (2.14) coincides with the high-pass filtered field in 
(2.2), so that —u'u' represents a 'fluctuation' contribution to the subscale stress, while 
q represents a 'systematic' contribution from the spatially-averaged, positive-definite, 
tensor product of velocity- increments. 



3. Convergent Expansion in Scale and Space 

In this section, we shall develop our convergent expansion for the turbulent stress. The 
key to this convergence is the locality of the stress both in scale and in space. Therefore, 
we shall discuss in turn these two properties, develop from them the resulting expansions, 
and establish their convergence properties. 

3.1. Locality in Scale 

Scale-locality is the property that only modes from length-scales near the filter scale I 
contribute predominantly to the stress. Recently, we have given a rigorous proof of this 
property, assuming only the inertial-range scaling laws that are observed in experiment 
and simulations (Eyink (2005)) and we refer to that work for a more complete discussion. 
Here we just recall that locality properties were demonstrated there by introducing a 
second 'test filter' F and an additional small length-scale 5 < I. A low-pass filtered 
velocity at scale S was then defined by 

u >s (x) = J d d r r s (r)u(x + r). (3.1) 

and likewise a stress contribution r >s , arising only from modes at length-scales > S, by 



t >s = u> 5 u> s -u> s u> s . (3.2) 
The property of ultraviolet (UV) locality of the stress is that 

lim t >s = t. (3.3) 

8^0 
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This limit means that modes at extreme subfilter scales (6 <C £) make little contribution 
to the stress at scale I. The result (3.3) was proved in Eyink (2005), with convergence in 
a strong L p -norm sense for any p ^ 1, under suitable spatial regularity assumptions on 
the velocity field. (In addition, some very mild moment-conditions must be satisfied by 
the filter kernels G and -T; see Eyink (2005).) A sufficient condition is that the scaling 
exponent £ 2 p of the (absolute) (2p)th-order moment of the velocity-increment 8u should 
satisfy the bound C, 2p > 0. There is also a similar property of 'infrared (IR) locality', 
which requires an oppposite condition ( 2p < 2p. Sec Eyink (2005). However, we shall not 
need to make use of IR locality in the present context. 

3.1.1. Multi-Scale Decomposition 

We can now reformulate the UV locality property as a multi-scale expansion of the 
stress tensor. First, we consider a corresponding multi-scale decomposition of the velocity 
field itself. Let us chose some parameter A > 1, e.g. A = 2 will be our standard choice. 
Then consider a geometric sequence of lengths £ n = \~ n £, for n = 0, 1, 2, .... with £o = £ 
and C\0asn^oo. For each of these we can define the corresponding low-pass filtered 
velocity field = it >£ ™, or 

u in) {x) = j d d r r ln {r)u{x + r). (3.4) 

This field includes modes at all length-scales down to £ n . We can also define a contribution 
to the velocity that arises, roughly speaking, from the length-scales between £ n -\ 
and £ n , by 

U M =«(")_ n>l. (3.5) 

It is convenient to set = u^K In that case, u^"' = J2k=o u for n ^ and, taking 
the limit as n — > oo, we get: 

oc 

u = J2 u[n] - ( 3 - 6 ) 

n=0 

This is the multi-scale decomposition of the velocity field. It is closely related to other 
similar scale decompositions, such as multiresolution expansions in wavelet bases, Paley- 
Littlewood decompositions, etc. Kraichnan (1974) used a multi-scale decomposition de- 
fined by banded Fourier scries in order to discuss locality properties of the turbulent 
cascade. The series (3.6) will converge in an L p norm, if, for example, the pth-order 
scaling exponent ( p of the absolute velocity-increment |£u(r)| is positive. In fact, in that 
case the series (3.6) has at least a geometric rate of convergence. 

The scale locality proved in Eyink (2005) can be restated in the present terms by 
defining the stress = -r >£ ™, or 



r (») = u (n) u (n) _ u (n) u (n) ) (3.7) 

which includes the contributions from all length-scales > £ n . If we substitute the ex- 
pansion = Y^k=o u ^ an d take the limit n — > oo, then we obtain a doubly-infinite 
series 



oo 



n=0 n'=0 

with t'"'" ' = uNuM — uW «["'!. This is the desired multi-scale expansion of the 
stress tensor. The term t'"'" ' represents a stress contribution from one velocity mode 
at length-scale £ n and another at length-scale £ n >. The series (3.8) converges absolutely 
in the L p -norm (and, in fact, at a geometric rate) under the same condition mentioned 
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in the UV locality statement, namely, the positivity of the scaling exponent of order 2p, 
(2 P > 0. This is a direct consequence of the concrete estimates in Eyink (2005). 

A remark should be made concerning the limit n — > oo in the expansions (3.6) and 
(3.8) above. For a finite (but arbitrarily large) Reynolds number, this limit need not be 
taken, practically speaking, because the stress contribution from the scales below the 
dissipative microscale is negligibly small. Instead, the series can be truncated at some 
n = rid corresponding to the length-scale of the viscous cutoff. However, for precisely 
this reason, there is also no difficulty in taking the limit n — > oo at a finite Reynolds 
number. In fact, the convergence rate of the expansion for n in the dissipation-range of 
scales is greater than in the inertial-range. As shown in Eyink (2005) the UV-locality of 
the stress depends upon the condition that the Holder exponent of the velocity satisfy 
a > 0, and, the larger a may be, the better this property holds. Since the velocity field 
is smooth in the dissipation range (a = 1), the UV-locality property is correspondingly 
stronger there than in the inertial-range where a < 1. 



3.1.2. Leading Terms and Strong UV-Locality 

Truncation of the series (3.8) at its leading term corresponds to making a strong UV- 
locality assumption. In that case, the exact stress is approximated as 

t w T m = uPltiH — vM vM. (3.9) 

This approximation neglects all interactions with modes at sub-filter lengths < £q = £. 
However, all interactions are retained between the modes at length-scales above the filter 
scale (both scale-local and IR scale-nonlocal ones). This approximation achieves closure 
for the stress in terms of ui ', essentially what is called the Similarity Model or Bardina 
Model in the LES literature (Mcnevcau & Katz (2000)). To simplify our discussion of 
this leading-order approximation, we may use a special notation for the low-pass filter 
of the velocity at scale I with respect to the test kernel r, namely, u = Fi * u. Thus, 
u and are different notations for the same quantity. We introduce also a simplified 
notation for the complementary high-pass filter, u" = u — u. 

The approximation (3.9) is rather extreme and the results in Eyink (2005) show that 
sub-filter scale modes can give a non-negligible contribution to the stress. Without re- 
peating all the details from that work, we would like to consider here briefly the mag- 
nitude of the error in making the strong UV-locality assumption. This approximation 
corresponds to replacing the exact velocity increment Su(r) in the formula (2.11) for the 
stress by Su^°\r) = Su(r). Let us assume that the velocity field has a Holder exponent 
a, that is, 8u(r) = 0(r a ) with < a < 1. (This notation means, as usual, that there 
exists a constant A so that |<5tt(r)| ^ Ar a , or, in a more dimensionally correct form, 
|<5it(r)| ^ CU(r/L) a , with L the integral length, U the rms velocity, and C a dimen- 
sionless constant.) Then it is not hard to prove that u"(x) = 0(£ a ); see Eyink (2005) 
for details. An error of this magnitude is made for all r in replacing 5u(r) by 5u(r). 
However, Su(r) = Su(r) + Su"(r) = 5u(r)[l + 0((£/r) a )} w Su(r) for r > I . Therefore, 
the substitution of 5u(r) for 5u(r) will be relatively accurate for increments over large 
separations r 3> £■ Of course, this substitution will not be accurate in the opposite case 
of small separations, because u is smooth and thus Su(r) w 0(r) -c Su(r) w 0(r a ) 
when r -C I. Thus, a relatively large underestimate results when the strong UV-locality 
assumption is applied to velocity-increments at small separations. 

Similar results hold for higher-order truncations, for example, for defined by (3.7) 
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with n ^ 1, or equivalently, 

n n 

- ( " ) = EE r[fe ' fc ' 1 - ( 3 - 10 ) 



k=0 k'=0 



This approximation assumes also UV-locality, but more weakly. It corresponds to replac- 
ing the exact velocity increment 5u(r) in the formula (2.11) for the stress by Su^ n \r), 
with an error 0(£"). This substitution will be relatively accurate when r > t n , but not 
for r < £ n . Here let us note that (3.10) does not yield a closed formula in quite the same 
sense as does (3.9), since it involves all of the components for k = 0, ...,n. If one 
assumes that G = T, then = u = u. However, even if G = r, one cannot in general 
obtain the higher terms «W, ...,u^ uniquely from knowledge of u. 

3.2. Locality in Space 

The turbulent stress is a priori non-local in space. The formula (2.11) expresses the stress 
as an average of velocity increments over separation vectors, which involves points, in 
principle, arbitrarily far away from the considered point. Nevertheless, the stress has some 
spatial locality properties, by virtue of the UV scale-locality discussed in the previous 
section. The latter property allows one to replace u by u^ n \ with an error 0(£") that 
becomes arbitrarily small for large enough n. The origin of the localness in space is 
then the smoothness of the filtered velocity field u^ n \ which will be even analytic if the 
filter kernel r has a compactly supported Fourier transform. This smoothness allows one 
to represent filtered increments, like 5u^ n \r), by a convergent Taylor-expansion in the 
separation vector r. The result is a formula that involves local gradients of the filtered 
velocities, i.e. velocity-gradients at the point where the stress is to be evaluated. 

3.2.1. Gradient Expansion 

Here wc develop the gradient expansion for the stress, once contributions have been 
omitted from arbitrarily small-scales. Let us first note the corresponding expansion of 
the filtered velocity itself. This field is smooth and thus the Taylor polynomial of 
degree m 

m pi Pd m . 

5u^ n - m \r-x) =J2 J2 ri \ ' " Td \ (di 1 • • • d p d d u^)(x) = -(r-V) p u^(x) 
p = 1 p 1+ ... +Pd=p P 1 ----P d - P=1 P- 

(3.11) 

converges to the increment 5u^ n \r;x) as m — > oo. If the Taylor polynomial <W"' m ) is 
substituted into (2.11), then it yields 



T (n,m) 



= J d d rG t {r)5u {n ' m) (r)5u {n ' m \r)- J d d r G e (r)5u {n < m \r) J d d r Gi(r)Su (n ' m) (r), 

(3.12) 

which is our basic approximation to the stress. 

An explicit expression is simplest to derive if one assumes that the filter kernel G is 
spherically-symmetric, as we shall do hereafter. In that case, averages over the directions 
of the increment vector r can be evaluated by a standard formula for averages of a product 
of an even number of vector components over the unit sphere in d space dimensions. The 
result, which is easily proved by induction, is that 

L-i W{dn)ni ^ • • • = d(d + 2)...[d + 2(p-l)] {i , .„ ,,, } *™ ' ' ' Si '^ 

(3.13) 
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where summation is over all of the {2p— 1)!! pairings {i' p , ip} of the 2p indices 

i\, 12, i2p-i, i2p- An average of a product of an odd number of unit-vector components 
is equal to zero. Particular cases are for p = 1: 



and p = 2 : 



/ w{dri)nirij = (l/d)6ij. (3-14) 



/ w{dn)ninjn k ni = 1 [5 i3 5 k i + SikSji + SuS jk }. (3.15) 
J S d-i a(d + Z) 



We thus obtain a stress approximation for m = 1 : 



T, 



lJ d dx k dxk 

where C2 = J d d r G(r) \r\ 2 is the 2nd-moment of the spherically-symmetric filter function 
G. Likewise, for m = 2 : 

r (n,2) = G2 e 2 du i ° U i + C 4 £i O Uj ° U j 



11 d dxk dxk 2d(d + 2) dxkdxi dxkdx, 



+ d ' C t*, ( ? + * C ® * A ^ A tiS^, (3.17) 



4d 2 (d + 2) 

where C2 is as before and C4 = J d d rG(r) \r\ 4 is the 4th-moment of the spherically- 
symmetric filter function G. In these expressions we may further substitute the multi- 
scale decomposition = X)fc=o u ^ to obtain expansions such as for m = 1 



[i'l 



^ - -J l^l^-Qx^-dx^ (3 ' 18) 

i=0 l'=0 



and similarly for t^ 71 '™^ with m > 1. Thus, we obtain a multi-scale gradient (MSG) 
expansion of the stress simultaneously in scale and in space. 

In Appendix A we prove that r("' m ) converges to in the limit as m — > 00. To 
keep the proof simple, we establish convergence in the spatial L^norm, requiring just 
finite energy for the velocity-field u. We assume also for the filter kernels that G decays 
faster than exponentially in space and that the Fourier transform r has compact support. 
These specific assumptions can doubtless be modified in various ways, but they simplify 
the details of the proof. From our discussion of scale-locality in the preceding section, we 
recall that also converges to r in the L 1 -norm as n — > 00, if the scaling exponent of 
the 2nd-order structure function satisfies (2 > 0. Thus, under these various assumptions, 
T (n,m) converges in the L 1 -norm to the exact stress r in the double limit taking first 
m — > 00 and then n — > 00. 

It is rather rare to be able to show that a systematic turbulence approximation scheme 
is convergent. For example, Kraichnan (1970) has discussed previous attempts to con- 
struct expansion schemes based upon Reynolds number, where convergence has proved 
quite elusive. A case in point is the gradient-expansion for the subscale stress proposed in 
Yeo & Bedford (1988), Leonard (1997) and Carati et al. (2001), based upon dcfiltcring. 
As we discuss in Appendix B, there are many fluid velocity fields with finite energy — even 
infinitely smooth ones — for which the expansion (B 4) does not converge. The problem 
becomes more severe the more rapidly the Fourier transform of the filter kernel, G(k), 
decays at large k. Our present study was motivated, in part, by the desire to overcome 
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this difficulty. As we have shown (Appendix A and Eyink (2005)), the multi-scale gradi- 
ent expansion that we have elaborated does indeed converge, under realistic and rather 
mild conditions on the turbulent velocity field and with very modest regularity assump- 
tions on the filter kernels. However, a price has been paid for this achievement. Unlike 
the expansion in Yeo & Bedford (1988), Leonard (1997) and Carati et al. (2001), our 
multi-scale gradient expansion is not closed in terms of the filtered field u — v}°\ but 
involves also the subscale fields u^\u^ 2 \ ... Thus, closure of our expansion requires an 
algorithm for estimating these unknown fields. 



3.2.2. Leading Terms and Strong Space- Locality 

Truncation of the gradient expansion r ( -™' m ' ) at small m values — e.g. approximating 
T («) ~ T ( n ^) to first-order in gradients, as in (3.18) — corresponds to making a strong 
space-locality assumption. If the expansion is truncated as well at small values of n, then 
both UV scale-locality and space-locality are assumed in a strong sense. E.g., setting 
n = 0, m = 1 gives 

which is the standard first-order Nonlinear Model for the stress (Leonard (1974), Leonard 
(1997), Borue & Orszag (1998), Meneveau & Katz (2000)). This observation gives some 
insight into the physical approximations underlying that model. 

We would like to make an estimate of the error involved in truncating the gradient 
expansion to a given order m of space-gradients. As in our discussion of scale-locality, we 
shall assume that the velocity field is Holder continuous, so that Su(r) = 0(r a ) for some 
< a < 1. Then it is not hard to show that the pth-order term in the Taylor expansion 
(3.11) of 8u^ m \r) scales as (r-V) p u^ = 0(£%(r/£ n ) p ) for each p > 1. For example, 
see Eyink (2005). Compared with the exact increment Su(r), we see that each term is 
an underestimate for r < £ n and an overestimate for r > £ n , and the error is greater for 
larger p. Truncated to a given small order m, the Taylor approximation has the correct 
order of magnitude only for r w £ n . If we sum over all values of p, to infinite order in m, 
then we recover the approximation Su^ (r), which we have seen is an underestimate for 
r < l n but relatively accurate for r > £ n . 

Although this gradient expansion converges, there is no small parameter involved (ex- 
cept for r <C £ n ). The series (3.11) converges only because of the inverse factorials l/p\ 
that make coefficients of higher-order terms quite small. Thus, we can expect that very 
large values of m will be required to make 8u^ n ' m \r) w Su^(r) w 5u(r) for r 3> £ n . 
Because the filter Ge(r) is assumed to decay very rapidly, increments with r » £ give 
little contribution to the stress and thus their poor approximation is not an issue. How- 
ever, increments for separations £ n < r < £ will give a significant contribution. In this 
interval the effective expansion parameter, r/£ n , takes values in the range from 1 for 
r = £ n up to £/£ n = A" for r = £. For small n, say n = or 1, the expansion parameter is 
always 0(1) in the relevant interval of r, and the series converges quite rapidly. However, 
as n increases, the rate of convergence in m degrades rather quickly. A crude estimate 
of the size of m required for an accurate approximation is m > A", in order for the 
terms 0(X np /p\) to be small for p w m. This estimate is probably too pessimistic, since 
it ignores cancellations that will occur in the average over r (cf. eq.(3.13)). However, we 
can be sure that the m required to obtain r^"' m ^ w will increase with increasing n. 
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4. The Coherent-Subregions Approximation 

We have proved that liirim^oo i-("' m ) = T^ n \ but also argued that larger orders of 
space-gradients m are required to achieve this limit for increasing n. This only stands 
to reason, because increasing the scale index n corresponds to adding finer small-scale 
structure to the velocity field. As the velocity field becomes rougher, higher-order terms in 
the spatial Taylor expansion become necessary in order to represent velocity-increments 
accurately across fixed separations. On the other hand, even a first-order expansion 
of increments, 5u^ n \r) rts {r-V)u^ n \ is correct on order of magnitude for separations 
r w £ n . Therefore, one should be able to get a reasonably accurate approximation by 
low-order gradients, if one uses such expansions only for this range of separations where 
they are order-of-magnitude correct. In the present section, we shall use this strategy 
in order to construct an approximate multi-scale representation rl"' m - ) for the turbulent 
stress. Although this modified expansion is no longer convergent to the exact stress, it 
may be more practically useful than the systematic approximation T ("> m ) ; because it 
should be more accurate for smaller orders m. 

Although it turns out not to be the most serviceable approach, a first natural idea is 
to represent increments r of length r > („ as the sum of end-to-end increments across 
separations of length l n and then to Taylor expand each of the individual increments. 
Thus, defining the unit vector r — r/r, one could write 

K K 

(5M (n) (r; x) « ^ Su^(£ n r; x + k£ n r) w £ n ^(f-V)u (n) (a; + k£ n r) (4.1) 

fc=0 k=0 

where K is the greatest integer less than or equal to r/£ n . The formula (4.1) is likely 
to be fairly accurate, since the expansion parameter is 0(1) for each term in the sum. 
However, this expression involves velocity-gradients evaluated at points x + k£ n r on 
spheres of radius k£ n about x for k = 0, 1, A" and spatially local expressions do not 
result for integrals over r when (4.1) is substituted into formula (2.11) for the stress. In 
order to get a simple, local expression, we should instead Taylor-expand always about 
point x. Thus, this approach does not lead to the desired result. However, at least it 
shows that an accurate representation of the stress is possible entirely in terms of filtered 
velocity-gradients of low-order, although the representation is spatially non-local. 

To obtain a local representation, we use Taylor expansions around x, but only for 
displacements where they are both rapidly convergent and accurate. Let us decompose 
the integrals in (2.11) into contributions from 'shells' 

S k = {r : 4-i > |r| > 4}, k > 1. (4.2) 

See Figure 1. Let us also define an 'outer shell' 

So = {r : |r| > 4} (4.3) 

and 'balls' 

B k = {r : |r| < 4-i}, k > 1. (4.4) 
formed from unions of the 'shells' Sk' with k' ^ k. From our earlier discussion we expect 
that, for r e Sk, 

Su(r) w <5«[ fc l(r) w 5u^^(r). (4.5) 

for any Taylor polynomial of degree m, since the expansion parameter is here r/4 ~ 1- 
Thus, we can obtain a rapidly convergent Taylor series expansion if we replace Su(r) by 
<5 M ( fe )( r ) f or r & Sk- However, this replacement implies that modes at length-scales < 4- 
are now represented only for the increments with r £ Bk- Furthermore, this ball Bk of 
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Figure 1. Shell Decomposition of Integration Region. The figure illustrates the integration region 
over increment vectors in the stress formula (2.11). Most of the contribution comes from the 
ball of radius £ — i around the point X. This ball is decomposed into 'shells' <Si, 1S2, S3, .... The 
Taylor expansion of the fc-scale contribution to the stress is accurate and rapidly convergent in 
the shell Sk- However, these shells have a smaller portion of the total volume as k increases. 

radius £k-i occupies only a fraction ~ \~ kd of the total volume of the region B\ which 
effectively contributes to the stress. Therefore, such a replacement omits important sub- 
scale contributions to the stress. To compensate for this, we can use the calculated stress 
contribution from the shell Sk to estimate crudely the missing part, by multiplying the 
calculated contribution with an enhancement factor of Nk = \ kd ■ This factor represents 
the number of subregions of volume <~ if inside the ball B\ of radius I. Multiplying each 
fc-scale contribution in the shell Sk by the factor Nk amounts to the assumption that 
each of the subregions gives a similar or 'coherent' contribution to the stress. We shall 
therefore call this heuristic estimate the Coherent- Subregions Approximation ( CSA ). Let 
us proceed to develop it more systematically. 

It is useful here to employ the formula (2.12) which decomposes the stress into 'system- 
atic' part g and 'fluctuation' part —u'u'. These can be represented further by multiscalc 
decompositions 

00 00 00 

fc=0fc'=0 fe=0 

analogous to (3. 6), (3. 8). Thus, g^ k ' k 1 represents the contribution from one velocity mode 
at length-scale Ik and another at length-scale ■ The multiscale decomposition of g can 
be re-organized as 

00 

q = E e m ( 4 - 7 ) 

fc=0 

where 

00 

g [k] = g [k,k] + J2{g^ + g^}. (4.8) 

i>k 

The term represents the contribution arising from a pair of modes, at least one at 
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length-scale Ik and the second at an equal or smaller length-scale. On the basis of these 
decompositions we can develop the desired approximation. 

First, let us consider the 'systematic' part g. The contributions g^- k,k ' scale as 0(£^£^i), 
when the velocity field has Holder exponent a. It follows that the dominant term in g^ 
is the first one on the righthand side of (4.8), or Qt k > k \ It is possible that the remaining 
terms sum to a contribution of similar order. However, the first term, £>[ fc ' fc l , is the average 
over space of positive-definite matrices, whereas the remaining terms have no definite sign 
and cancellations can be expected in the summation in (4.8). Thus, we expect that 

g [k] M g [k,k] = J d d rGl ( r ) du W(r)6uW(r). (4.9) 
This motivates us to define the CSA value of g^ k \ mth-order in gradients, as 

g lkUm) =N f d d rG ^ r ) §u [kUm)^ §u lk],(m)^ (41Q) 

As discussed earlier, the factor on the righthand side of (4.10) corresponds to making 
the assumption that each subregion of B\ gives a contribution to the integral (4.9) similar 
to that of the fcth 'shell' Sk- This is reasonable, since the integrand is positive-definite 
and thus there will be little cancellation between the contributions from the different 
subregions, which can be expected to add together coherently. Using then the replacement 
(4.5) in the integral over Sk gives (4.10). As an additional argument in favor of the 
enhancement by Nk, let us note that g^^" 1 ^ defined in (4.10) with this factor gives an 
0(t\ a ) contribution to the stress, of the correct order of magnitude. See Appendix C. 

Similar considerations apply also to the 'fluctuation' velocity it'. The /cth-scale contri- 
bution u' [ fe l in (4.6) is written exactly as 

w'[ fe l =- J d d rGi(r)Su^(r). (4.11) 
We then propose the CSA value of u'^ k \ mth-order in gradients, as 

^[kUm) = _ N l/2 f d d rGl{r)5u [k],(ra) {r) _ (412) 
JS k 

1 /2 

Note the change in the enhancement factor to N k . The integrand in (4.11) has no 
definite sign and the Nk different subregions should not be expected to contribute co- 
herently. In a work on analytical closures, Kraichnan (1971) made a similar argument 
about the shear contribution from small scales, writing that 'random cancellation effects 
over the domain 1/k in linear dimension should reduce the effective shear of the high 
wave- numbers according to the ^/~N law.' Analogous reasoning motivates us to multiply 
the righthand side of (4.12) by N^ 2 . 

This set of approximations altogether yields the CSA-MSG expansion for the stress, 
nth-order in scale index and mth-order in gradients: 

n n 

T (n,m) = ^ C W.("0_ u} kUm) u} k ' Um) (4.13) 

fe=0 fe,fe'=0 

Using the results for m = 2 as illustration, we can write 



d dxi dxi 2d(d + 2) dxidx m dxidx v 
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(4.14) 



and 



;[fe],(2) _ 




(4.15) 



As in (3. 16), (3. 17), we have used (3.13) in order to average over the directions of sepa- 
ration vectors r in (4. 10), (4. 12). Note that (4.14) has the same form as (3.17) for t (m) 

except that the coefficients are different. The constants C p for p = 2, 4, ... are the partial 
pth-momcnts of the kernel G over the fcth 'shell' S k , multiplied by the factor \( d +P) k . 
Expressions are given for these constants in Appendix C, with G a Gaussian filter. 

Our rather rough estimates should obviously be taken with a large grain of salt and 
are intended to be accurate qualitatively, but not more than order-of-magnitude accurate 
quantitatively There is clearly ample room to improve the accuracy of the scheme, and 
many variants and refinements might be fruitfully considered. The basic approximation in 
the 'coherent-subregions' assumption, i.e. estimating missing small-scale contributions to 
the stress from their effects in subvolumes, will tend to enhance the level of fluctuations. 
However, the CSA stress rl™'™' in (4.13) is still likely to be superior to the systematic 
MSG expansion r^™' 1 ™-' when n > 1 and m is relatively small. The approximate stress 
T (n,m) conver g es ra pidly i n the limit m — > oo, to some value rl™' 1 which is, hopefully, a 
reasonable approximation of r^ n \ requiring only moderately large values of m uniformly 
in n. It achieves our goal of providing a local expression for the stress which involves 
only filtered velocity-gradients of low-order. However, like T^ n,m \ it is not a proper 
constitutive relation, because it is not closed in terms of u = vffl. Of course, it already 
makes testable predictions for the stress, if the smaller-scale velocity fields for k ^ 1 
are determined from experiment or DNS and then substituted into the model. We report 
results of such a study elsewhere. Nevertheless, an a priori closure procedure would be 
useful for modelling purposes. This could be accomplished by a stochastic mapping which 
estimated the velocity gradients Vm'^ , VVii' fc ' , etc. for k > 1 from the corresponding 
gradients for k = 0. We hope to make this the subject of a future work. 

5. The Multi-Scale Gradient Expansion in 3D 

As an application of the general scheme, we shall consider here the turbulent cascades 
of energy and helicity in three space dimensions. In a following work (Eyink (submit- 
ted)), we discuss the MSG expansion of the stress for the inverse energy cascade in two 
space dimensions. Many other applications can be considered, such as turbulent vorticity 
transport in the 2D cnstrophy cascade (Eyink (2001), Chen et al. (2003)), or the tur- 
bulent stress tensor and electromotive force in 3D magnetohydrodynamic cascades. The 
technical aspects of the expansion are similar in all of these cases. As we shall see in this 
section, our method yields a number of interesting predictions for the turbulent stress in 
3D that may be tested either numerically or experimentally. 



We shall confine ourselves here to considering just the first-order (m = 1) term in the 
gradient-expansion, or t^ 1 ^ in (3.16). We have already noted that this first-order ap- 
proximation is unlikely to be very accurate for larger n. On the other hand, because 
of the scale-locality of the energy cascade, only relatively small values of n need to be 
considered and thus a first-order approximation may be adequate. Furthermore, except 
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for the coefficient, t^™' 1 ' has the same form as the term g!™''^ in (4.13) for m = 1 : 

and, in addition, the 'fluctuation' term in (4.13) vanishes for m — 1. Thus, the first-order 
CSA expansion has the closely similar form 



(n,l) 



_ 1 s^H k h* dU dU (*>2) 

~~ d h 2 k ~teT~tei" { } 

We expect that t!™' 1 ^ in (5.2) for large n will be reasonably accurate in the 3D inertial- 
range. For example, the estimates in Appendix C show that the fcth term in (5.2) scales 
<~ 0(£l a ) when the velocity field has Holder exponent < a < 1. This is the correct 
order of magnitude for the contribution to the stress from scale k and illustrates the UV 
locality of the stress. The series in (5.2) then converges at a geometric rate in the limit 
n — ► oo and has the correct overall mag nitude ~ 0(£ 2a ). By contrast, in (3.16) 

does not have the correct order of magnitude as n — > oo, but is too large by a factor of 
(£/£ n ) 2 - This is due to an overestimate in r^™' 1 ) of velocity-increments at large spatial 
separations, arising from the first-order Taylor expansion. Thus, our concrete results 
below for t^™' 1 ) in 3D should be more properly reinterpreted, when n is large, for the 
approximation ri™' 1 ^ in (5.2) instead. 

In any case, we have in 3D the formula for the filtered velocity-gradient 

,(«) 

b ij ~ o e ijkU k (0.6) 



du i ' = _ £ , ,(n) 

dxj ij 2 



ij 3 2 



in terms of the filtered strain tensor S^ n \ the filtered vorticity vector a/") and the 
antisymmetric Levi-Civita tensor e^fe. If (5.3) is substituted into (3.16), it yields 

The separate terms in this expression have interesting physical interpretations. The first 
term is proportional to the strain-matrix squared: 

[S H ]2= ^| CT (n)|2 e (n) e (n) j (5 . 5) 
p=l 

where and are the eigenvalues and eigenvectors of the strain matrix 
satisfying a[ n ^ + + = 0. This term represents a tensile stress of magnitude 
(l/3)C2(<7p"^) 2 exerted along each of the principal strain directions for p = 1,2, 3. 
The last term in (5.4) quadratic in the vorticity likewise represents a tensile stress along 
the two directions orthogonal to the filtered vorticity. However, the first part of that term 
proportional to the Kronecker delta function is an isotropic stress or turbulent pressure, 
which does not contribute to the deviatoric stress. The other half of the term, propor- 
tional to bj^oj^ l \ is equivalent to a contractile stress of magnitude — (l/12)C2(w^™^) 2 
exerted along vortex-lines. Thus, one of the important effects of subscale modes is an 
induced tendency for lines of filtered vorticity u/ n ) to resist lengthening. This 'elastic 
response' of vortex-lines is well-known in other contexts — for example, turbulence under 
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rapid-distortion (Crow (1968)). However, the most novel of the stress terms in (5.4) is the 
middle one, which is given by a certain 'cross product' of strain and vorticity. More pre- 
cisely, we have defined (u/ n ) xS (ri) )„- = e^wj^S^ . Note that this tensor is orthogonal 

to the strain at the same scale, S'"':(ti;'"'x5'' 1 ') = 0, so that we call it the skew-strain. 
The middle term of (5.4), proportional to this skew-strain, is a sum of shear stresses 

1 3 

-C 2 ]T U WaW sin0(") [e^e^ + e^e^] (5.6) 

P =i 

where 9^ is the angle between u/") and e p n ^ and e p n ^ is the unit vector orthogonal to 
both u/™) and e p n \ given by the righthand rule. If we introduce the new unit vectors 
e p± = I e p"' ) ^ ep™^]/\/2, then (5.6) becomes 

l(7 a £y»>aW sin^)[e^e^ - e p n }e p n }\ (5.7) 

Hence, there are both tensile and contractile stresses exerted along the vectors ej^. 

These are obtained from the strain eigenvector e p n ^ by rotating it ±7r/4 radians around 
the normal component of the vorticity vector a/™'. 

The above vector formalism helps to make clear the geometry of the various stress 
contributions. However, it is perhaps more conventional to write these stresses in terms 
of the fluid deformation matrix D^ n \ defined by — duf 1 ^ /dxj , and its symmetric 
part 5^ and anti-symmetric part f2^ n \ Of course, S*-"-* is the strain tensor and l?*-™' is 
related to the vorticity vector u/ n ) by the standard relation Q^ 1 " 1 = — {1 / 2)e i j k uj^ L \ In 
terms of the deformation matrix the lst-order term (3.16) in the MSG expansion can be 
written (in fact, in any dimension d) as 

T< n,1 > = ^C2l 2 D (n) [D in) } T - (5.8) 
The decomposition analogous to (5.4) is then 

T (n,i) = \ C2 p |s(«) S («) + [^(™),S (n) ] - (2 {n) Q {n) } (5.9) 

where [Q {n \ S (n) ] = Q {n) S {n) -S (n) a (n) is the commutator matrix. Thus, in 3D, S {n) S {n) 
is the strain-squared as in (5.5), -/?("> fl< n > gives the tensile stress in the plane normal 
to vortex lines, and [Q^ n \S^ n>> ] is the 'skew-strain.' 

5.2. Energy Cascade in 3D 

It is interesting to consider the consequences of the stress in formula (5.4) for the energy 
cascade. When (5.4) is substituted into equation (2.8) for the energy flux, one gets the 
following result to first-order in gradients: 

i7("' J ) = |c 2 ^|-Tr (S(S (K) ) 2 ) +^-(w(")) T S(a;(")) + S:(S ( " ) Xu;("))|. (5.10) 

The middle term has an obvious physical meaning. It represents the rate of work done 
by the filtered strain S in order to stretch the lines of vorticity a/™) against the resisting 
contractile stress of the subscales. This remarkable relationship between energy flux and 
vortex- stretching was already observed by Borue & Orszag (1998) using the Nonlinear 
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Model, which is a special case of our result for n — and G = T: 

il(°' i ) = ^ 2 {-Tr(sV^ T S^} (5.H) 

Vortex-stretching was suggested long ago by Taylor (1938) as the basic dissipation mech- 
anism of 3D turbulence. Note that the first term in (5.11) proportional to the strain 
skewness can also be related to vortex-stretching, on average, using a relation of Betchov 
(1956). His result states that for an incompressible fluid — (Tr (5 )) = (3/4)(tJ T 5uJ), 
where (•) denotes either ensemble average over a statistically homogeneous turbulence 
or any volume-average where boundary terms from integration-by-parts can be ignored. 
According to Betchov's relation, precisely 75% of the mean energy flux in the Nonlinear 
Model comes from strain skewness and 25% from vortex-stretching. It is interesting that 
Betchov's relation can be generalized as follows: 

-(Tr (s(S (n) ) 2 )) = i((aJ) T S("V n) )> + \{(^) T S(u:^)). (5.12) 

This result is proved in our Appendix D under the same assumptions as Betchov's. By 
using (5.12) the first term in (5.10) can be related to vortex-stretching in general for 
all n. However, the last term in (5.10) appears to be fundamentally different. It appears 
only due to contributions of subscales and there is no analogue in the Nonlinear Model 
(5.11) for n — 0. Some additional insight on that term can be obtained by rewriting 
S:(S (n) Xw (n) ) = u>(")-(SxS^), where SxS (rl) is the dual vector corresponding to the 
antisymmetric commutator matrix [S, S^], i.e. (SxS^)j = (l/2)ejjfe([S, S^])jfc. Thus, 
this new term arises from rotation of the subscale strain relative to the filtered strain 
5, and vanishes if the orthogonal cigenframes of these two symmetric matrices coincide. 
Using the CSA stress in (5.2) one gets a similar result as (5.10): 

^ = | fy^ 1 " {~ Tr (^ SW ) 2 ) + ^ [fcl ) T 5(^ [fcl ) + S:(S[ fc lxu;[ fc ])| . (5.13) 

The remarks we have made above on physical interpretation of Tjf™' 1 ) apply equally here. 
However, the sum in (5.13) also has a limit for large n and we expect that it gives a quite 
reasonable model for energy flux in 3D. Using the generalized Betchov relation (5.12) the 
CSA mean flux can be written as 

<tfi B,,) > = {'-{{wys^)) + |<(«W) T 3(«W)> + @(SMxa,M)>} . 

k = ^ ' 

(5-14) 

The first two terms in each summand arise from vortex-stretching and will tend to be 
positive, certainly for small k when S' fc ' cx S, u:^ oc uJ. On the other hand, the third 
term from skew-strain then nearly vanishes. For the latter term to be important, there 
must be some characteristic rotation of S' fe ' relative to 5 as k increases. Of course, it is 
not hard to see that each term is zero on average when S^,u;t fc l are uncorrelated with 
5, which must be expected in the limit as k — > oo. Therefore, it only for intermediate 
values of k that the third term can contribute to mean energy flux. 

The physical mechanism of energy cascade by these stress terms can be illustrated by 
the following: 

Example 1. Vortex Tube Stretched by a Constant Strain 

We consider a small-scale cylindrical vortex-tube, parallel to the z-axis, with circular 
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Figure 2. Energy Cascade by Vortex-Stretching. The figure illustrates a cylindrical tube of 
parallel vortex lines (center) in a constant strain field with stretching direction along the vortex 
axis (left). The result, shown on the right, is that the vortex is stretched and its cross-sectional 
area shrunken. The 'spin-up' of the vortex increases its energy and generates more positive 
(tensile) stress in the plane perpendicular to the vortex axis. 



cross-section of radius R — l k and with vorticity magnitude cj = u^. This is an exact 
stationary solution of the 3D Euler equation with two-dimensional symmetry. Thus, it 
may be described by a pseudoscalar stream function 



^ k \x,y) 



(l/4Vo[i? 2 -r 2 ] 
-(l/2)u R 2 ln(r/R) 



r < R 
r>R 



(5.15) 



where r is the radial distance from the z-axis in cylindrical coordinates (r,8,z). The 
corresponding velocity field is 



u 



Ik] 



(x,y) 



(l/2)u> zxr 
(l/2)Lu (R/r) 2 zXr 



r < R 
r > R 



(5.16) 



where z is the unit vector in the z-direction. See Figure 2. This small-scale field is now 
superimposed with a large-scale velocity 



u 



-{a/2)x 
az 



with deformation matrix D 



(5.17) 



(5.18) 



This is a pure large-scale strain D = S with vorticity ZJ = o. If a > 0, then this cor- 
responds to an axisymmetric stretching along the z-direction and compression in the 
other two directions (Figure 2). The combination of the large-scale and small-scale fields 
gives an exact solution of the 3D Euler equation (dt + u-V)u>^ (t) = (u>^(t)-'V)u = 0, 
where u = u + u^^t) and where u>^(t) is the same as the initial vorticity field, made 





' -a/2 








D = 





-a/2 













a 
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time-dependent by the substitutions uj(t) = e CTt o>o and R(t) = e~ at l 2 R. For example, 
see Neu (1984), who considers a more general set of solutions. In the present case, the 
small-scale vortex is stretched along its axis and, by incompressibility, its cross-section 
shrinks. To conserve the circulation around the tube, the vorticity and the velocity in 
the small-scales correspondingly increase. See Figure 2. This 'spin-up' by the large-scale 
strain results in a transfer of energy to the small-scales. 

The process can be understood from our general formulas above. Without loss of 
generality, we can focus on the instantaneous transfer at the initial time t = 0. The 
velocity-gradient tensor in the small-scales is then 



D [k \x,y) 



,[k] 



,[k] 



-1 

1 


sin(20) 
- cos(26>) 








- cos(26») 

- sin(26») 





r < Ik 



r > 



(5.19) 



This is purely rotational for r < l k and is a pure strain for r > l k . Substituting into (5.1) 
gives the stress 



7* 



1 
1 




1 

(4/r) 4 



r < 4 
r > 4 



(5.20) 



to first order in gradients. This result represents the — term for r < £ k and the 

term for r > l k . The 'skew-strain' vanishes identically for this right cylindrical 
vortex tube. We see that the net stress is tensile in the 2D plane perpendicular to the 
vortex tube, set up by the velocity circulating around the vortex axis. Its deviatoric 
part includes a contractile stress along the vortex axis. These stresses oppose the axial 
stretching and lateral compression by the large-scale strain, and increase in magnitude 
as the vortex spins up. The work of the large-scale strain against these resistive stresses 
is the basic mechanism of energy transfer to the small-scales. 

It is interesting to observe that if a < 0, then the energy flux corresponding to (5.20) 



n [fc],(D 



1 r<4 

(4/r) 4 r>4 



(5.21) 



is negative and the large-scale strain 'spins down' the small-scale vortex, by the time- 
reverse of the process considered above. What is crucial for foward energy transfer is 
that the vortex should align with a stretching direction of the large-scale strain rather 
than with a shrinking direction. We know from the relation of Betchov (1956) that, in an 
incompressible flow, mean vortex-stretching requires that there be typically two positive 
strain eigenvalues and one negative eigenvalue. This tendency has been confirmed for 
dissipation-range velocity gradients by DNS (Ashurst et al.(1987)) and for inertial-range 
(filtered) velocity-gradients by experiment (Tao, Katz & Meneveau (2002), Van der Bos 
et al. (2002)). Furthermore, these empirical studies have shown that the vorticity vector 
tends to align with the intermediate, weakly stretching eigendirection of the strain at the 
same scale, rather than with the strongest stretching direction. Some theoretical under- 
standing how this occurs can be obtained from simple Lagrangian dynamical models of 
the velocity-gradients (Vieillefosse (1982), Vieillefosse (1984), Cantwell (1992), Chertkov, 
Pumir & Shraiman (1999)). Thus, the situation in turbulence is slightly different from 
that which we imagined in our simple example above. However, the mechanism of the 
energy transfer process appears to be essentially the same. 
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5.3. Helicity Cascade in 3D 

It is well-known that 3D smooth solutions of the inviscid, incompressible fluid equations 
have in addition to the energy a second quadratic invariant, the helicity (Moreau (1961), 
Moffatt (1969)): 

H(t) = J d 3 x u(x,i)-u(x,t). (5.22) 

When helicity is input at large scales together with energy, then there is in 3D a joint cas- 
cade of both invariants to high-wavenumber (Brissaud et al. (1973), Kraichnan (1973)). 
The flux of helicity can be expressed quite similarly to the flux of energy in (2.7), as 

A = -2VlJ:t (5.23) 

See Chen, Chen & Eyink (2003). Formula (5.23) is quite intriguing, since it implies that 
the stress r must be correlated simultaneously with both the velocity-gradient and the 
vorticity-gradicnt in a joint cascade of energy and helicity. Our work sheds some light on 
how this is achieved. 

A vorticity-gradient may be decomposed into symmetric and anti-symmetric parts, as: 

Rij H~ ^ij Rij ~2^ijk^,k (5.24) 

where R= (l/2)[(Vw) + (V«) T ], 5 = (l/2)[(Vw) - (Vw) T ] and £ = Vxw. We may 
write the helicity flux also as A = —2R: r, because of the symmetry of the stress tensor. 
Therefore, if we substitute the first-order stress formula (5.4), then we get the expression 
for helicity flux analogous to (5.10): 

/1("< J ) = Jc 2 £ 2 j-Tr (7?(S (n) ) 2 ) +^-(a;(™)) T R(a;(™)) + R:(S ( " ) Xa;("))| (5.25) 

This is the exact expression to first-order in gradients. Of course, we can also write down 
a CSA expansion for helicity flux, 

A (n,i) = I fyfil j-Tr (R(SW) 2 ) + + R:(S [k] xJ k ^ , (5.26) 

analogous to (5.13) for energy flux. All the stress components — from strain-squared, 
from vortex contraction, and from skew-strain — contribute to the helicity flux. Note 
that the generalized Betchov relation from Appendix D can be applied to give 

-(Tr (7?(S [fc1 ) 2 )) = l(f£ k] u [k] ) + |((w w ) T & [tl ), (5.27) 

where (•) denotes a homogeneous average. Thus, the strain-squared contribution can be 
replaced on average with the above two terms. 

It is known that the helicity cascade is local in scale (Eyink (2005)). Therefore, it is 
interesting to consider the n = contribution, which, assuming G = J 1 , coincides with 
the helicity flux for the Nonlinear Model of the stress. Now, it is not hard to see that 

Z3 T RE3 = V- [(1/2)|ZU| 2 ZU] , (5.28) 

which is a total derivative. Thus, 

{l/A)(TD T Ru) = 0, (5.29) 
and the contractile stress along vortex lines gives no contribution to the UV-local part 
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of mean hclicity flux. Combining (5.29) and the generalized Betchov relation (5.27) gives 
also 

-(Tr{R(Sf)) = ±(fSu>). (5.30) 
Therefore, the total UV-local (n = 0) contribution to mean helicity flux is 

(A (0 ^) = - 3 C 2 l 2 {^(fSuJ) + (tt(SxuO)} (5.31) 

Equivalently, this is the Nonlinear Model expression for mean helicity flux. The first term 
arises from the stress proportional to strain-squared and the second term from the stress 
proportional to skew-strain. 

The two terms in (5.31) can be related by the following identity 

il T SuJ = -E:(Sxuf). (5.32) 

This is easily proved by substituting on the right Sij = — (l/2)ejjfc£ fe , then using the 
definition of the skew-strain and the identity e m ijtmki = SikSji — SuSjk- Since (Vu;) T = 
~R—~E, using (5.32) in (5.31) gives 

= |c 2 ^((VuJ) T :(SxuJ)) (5.33) 

Thus we see that both the symmetric and anti-symmetric parts of the vorticity-gradient 
can contribute to helicity flux. The result (5.33) also makes clear the important role of 
the skew-strain in the 3D helicity cascade. It is noteworthy that skew-strain makes no 
UV-local contribution to energy flux at all and is thus free to adjust as necessary to 
maintain the helicity flux in a joint cascade of both invariants. 

If we use the result (l/2)£ T S [fe W fc ] = - S:(S [k] xw^) analogous to (5.32), the gener- 
alized Betchov relation (5.27), and the expression (5.26), then we obtain 

<^)> = l±cftl {<(Vu7)T:(SMxu,M)> + I^M^M))} , (5.34) 

for the CSA expansion of mean helicity flux. This is analogous to the similar result 
(5.14) for the mean energy flux. However, note that it is now the first term which makes 
a UV-local contribution while the second only contributes for intermediate values of k. 

We now consider a simple example to illustrate the mechanism of hclicity cascade by 
these stress terms: 

Example 2. Vortex Tube Twisted by a Constant Screw 

We take as our model of the small-scales an exact stationary solution of 3D Eulcr 
equations which was previously considered by Moffatt (1969) as an example of a helical 
flow with continuous vorticity distribution. It is a two-dimensional but three-component 
(2D-3C) velocity field, closely related to our previous Example 1. Indeed, the horizontal 
components (u^ k \v^) of the velocity field are the same as those in (5.16), obtained from 
the 2D stream-function ip^ in (5.15). However, this is now supplemented with a vertical 
velocity component 

(1/4)ujop[R 2 -r 2 } r<R 
r>R 

It can easily be shown that the resulting total velocity field is a stationary Euler solution 
(e.g. see Moffatt (1969), Section 6(a).) If p = 0, then this solution coincides with that in 



wW(x,y) = \ WWl" - < J < (5 35) 
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Figure 3. Helicity Cascade by Vortex-Twisting. The figure illustrates a cylindrical tube of 
parallel vortex lines (center) in a constant screw field with twisting direction along the vortex 
axis (left). The result, shown on the right, is that the vortex lines are twisted into helices. 
The tilting of the small-scale vorticity vector and the solenoidal generation of an axial velocity 
create helicity at small-scales. A positive (tensile) stress occurs, in the plane perpendicular to 
the vortex axis at loose winding and along the twist axis at tight winding. 



our Example 1. The meaning of the parameter p can best be understood by considering 
the associated vorticity vector u)^ = (a^ , 0& , 7^). Of course, the vertical component 
7^] is the same as in Example 1, = ujq for r < R and = for r > R. The horizontal 
components are obtained using the vertical velocity as a 'stream function': 



dw^/dy 
-dwW/dx 



-y 



(5.36) 



for r < R and = for r > R. The vortex lines are helices winding around the z-axis with 
'pitch' Air/p, i.e. making one counterclockwise revolution in that vertical distance. It is 
not hard to check that the solution given by (5. 15), (5. 16), (5.35), (5.36) has a constant 
helicity density = u^-u>^ = (l/A)u^R 2 p for r < R and = for r > R. 
We shall take as our model of the large-scales the velocity 



u 



-pyz/2 
pxz/2 




(5.37) 



This corresponds to a constant screw, i.e. solid-body rotation in each plane parallel to 
the xy-plane with an angular velocity ~pz/2 that grows linearly in z. It is a righthand 
screw for p > and a lefthand screw for p < 0. This velocity field has both non- vanishing 
strain and vorticity: 



(5.38) 











-py/4 ' 




" -px/2 ' 


s = 








px/A 


, w = 


-py/2 




. -py/4 


px/4 







p z 
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Furthermore, the vorticity-gradient matrix is constant and symmetric: 



-p/2 







-p/2 




R. 



(5.39) 



It is not hard to check that the small-scale velocity field defined previously becomes 
an exact time-dependent solution of the 'rapid-distortion equation' (<9 t + it-V)u;M (i) = 
(o;[ fc l (£). V)w = 0, with w = u + u^ (i), if the parameter p = pt. (This is no longer equiv- 
alent to the full 3D Euler equation, since the large-scales have also vorticity, neglected 
here.) The small scales begin as the undisturbed vortex tube of Example 1, whose fila- 
ments are then twisted into helices by the large-scale screw. The resulting vorticity field 
of coiled helices generates the axial velocity of (5.35) by solenoidal action, producing 
a net helicity in the small-scales. See Figure 3. 

Let us now consider the helicity transfer process, based upon our general formulas for 
the stress. We shall only consider the space region r < R, since the small-scale fields 
outside the tube are the same as for Example 1. Inside the small-scale vortex tube the 
velocity-gradient tensor is 



0-10 
1 

— px —py 



(5.40) 



Substituting into (5.1) gives the stress 



-c lfc] M fe J 

12 1 



4| 2 



1 py 

1 — px 
py —px p 2 r 2 



(5.41) 



to first order in gradients. Unlike Example 1, all three terms in the stress formula (5.4) 
[or (5.9)] are present, including that from 'skew-strain'. When p is small (p£ k <C 1) then 
this is essentially the same result as in Example 1, but when p is large (p£ k ^ 1) the 
dominant stress is tensile along the screw axis. The resulting helicity flux has the form 



A [fc],(i) = lc [k] 
6 



p\J k h k \* [i-p 2 



(5.42) 



When p is small, there is a net transfer of helicity to the small-scales of the same sign 
as the large-scale screw. This arises from the weakly local transfer produced by the 
large-scale vorticity-gradient acting against the contractile stress along the small-scale 
vortex. However, for large p the sign of helicity flux reverses, as the more tightly wound 
vortex lines produce a net tensile stress along the screw axis by solenoidal action. The 
contributions to helicity flux of the separate S (n) S (n) , [Q {n) , S {r % and -/?<") I)W stress 
terms in (5.9) are 



±c [k] - P \pJ k k k \\\ -\c [k] - P \pJ k k k \\\ \c [k] 



P|w [fc] 4| 



1 - 



1 



2 2 

p r 



(5.43) 



respectively. We see that the flux comes mainly from the —Q^Q^ term for pi k -C 1, 
as claimed above. For pl k » 1 all three terms contribute, with the flux contribution from 
the 'skew-strain' six times bigger than that of either other term. 

This example also illustrates the cascade of energy. Indeed, as was observed already 
by Brissaud et al. (1973), the transfer of helicity necessarily involves also the transfer of 
energy. The formula (5.41) for the small-scale stress and (5.38) for the large-scale strain 
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yield the following result for energy flux: 

n [ fe ],(D = l_c lk] p P \J k k k \^ (5.44) 

There is a net forward transfer of energy if the signs of p and p are the same. In that case, 
the large-scale screw winds the small-scale helical vortex-lines more tightly, and kinetic 
energy is generated in the axial velocity component by the resulting solenoidal action. 
The contributions to energy flux of the separate S {n) S {n \ [Q {n \ S (n) ], and -Q {n) Q {n) 
stress terms in (5.9) are 

0, lc [fe Wk W 4|V, A^WM^IV, (5.45) 

respectively. Inside the vortex, there is no energy transfer from the strain-squared and 
instead equal amounts of the energy flux are due to the contractile stress along vortex 
lines and the stress proportional to 'skew-strain'. 



6. Conclusions 

In this paper we have developed a novel approximation for turbulent stress, via a multi- 
scale gradient (MSG) expansion. This scheme represents the stress by an expansion in 
scales of motion and in orders of space-gradients. A major result (Appendix A) is that 
this expansion converges and furthermore at a rapid rate for the 'strongly local' part of 
the stress from the resolved scales and adjacent subscales. However, the convergence of 
the spatial Taylor expansion is expected to be much slower for stress contributions from 
scales further below the filtering scale. Therefore, we have developed a more approximate 
expansion, which should give a reasonable result at all scales with just a few low-order 
velocity-gradients. This 'coherent-subregions approximation' (CSA) is based on the as- 
sumption that the velocity increments across all the subscale separation vectors should 
give a similar result, at a given scale, as those for separation vectors in a 'shell' where 
the gradient-expansion is accurate and rapidly convergent. 

An important application of our methods has been presented to the three-dimensional 
turbulent cascades of energy and helicity. Our main results are the formulas (5.4) for 
the stress, (5.10) for the energy flux, and (5.25) for the helicity flux, exact to first-order 
in gradients. We have also developed the corresponding CSA expressions, (5.2), (5.13), 
and (5.26), which are more heuristic but which should give a better representation of 
the very small-scale contributions than the exact first-order results. We have generalized 
Betchov's well-known relation, which relates mean vortex-stretching and strain-skewness 
at the same scale, to a similar relation between different scales ((5.12) and Appendix D). 
This relation allowed us to derive expression (5.14) for mean energy flux and (5.34) for 
mean helicity flux, and to analyze the expected contributions at different scales. Finally, 
we have discussed the physical mechanisms of energy and helicity cascade, in terms of 
our analytical formulas. We have shown by means of simple exact solutions of 3D Eulcr 
equations that our results are consistent with energy transfer by Taylor's mechanism of 
'vortex-stretching' and with helicity transfer by a mechanism of 'vortex-twisting'. 

There are many implications of the present work for experiment and simulation, for 
theory, and for modelling. 

A host of testable predictions have been provided for laboratory experiment and for 
numerical simulation by our detailed formulas for turbulent stress, energy flux, and he- 
licity flux. The expansions in scale and in space have been proved to converge, but the 
rate of convergence could be even faster than what has been rigorously established and 
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empirical studies can determine this. The very distant subscales have been proved in 
Eyink (2005) to give decreasing contributions to stress, and analytical closures make fur- 
ther quantitative predictions about the mean amount of energy and helicity transfer from 
each scale of motion (e.g. see Kraichnan (1971), Andre & Lesieur (1977)). Our multiscale 
formalism provides a convenient framework within which to check these predictions, par- 
ticularly for experimentalists who cannot easily calculate spectral transfers. As to the 
gradient-expansion, our first-order expressions (5.4), (5.10), and (5.25) should give good 
results for the strongly local contributions, without further approximations. Experiment 
and simulation can also check the validity of our CSA expansion, which is based upon 
a bolder assumption. Assuming that our results are empirically confirmed, experiment 
and simulation can also determine the relative magnitudes of the various terms in our 
formulas. This will help to shed further light on the detailed physical mechanisms which 
underlie the turbulent cascades. 

Our results suggest several further fruitful directions for theory. Previous dynamical 
models of velocity-gradients (Vieillefosse (1982), Vieillcfosse (1984), Cantwell (1992), 
Chertkov, Pumir & Shraiman (1999)) have investigated only alignments between objects 
at the same scale. However, as our results should make evident, alignments between 
velocity-gradients at distinct scales are also of great importance in supporting turbulent 
cascades and these inter-scale relations have received scant attention so far. Improvement 
of our various approximations is another important avenue for theory, particularly the 
CSA scheme, where many possible refinements are apparent. Finally, theoretical methods 
to estimate subscale velocity gradients from the resolved ones are strongly motivated 
by our work. This would lead to LES modeling schemes, similar to those reviewed in 
Domaradzki & Adams (2002), but based upon a clearer picture of the physical processes 
involved. If there is a universal mechanism which generates small-scales from large-scales, 
then modeling this generation process should lead to the most robust and generally 
applicable LES models. 

A virtue of our approach is its wide range of potential applications. Since it is based 
upon a very general feature of turbulent cascades — i.e. their locality in scale and in 
space — the same scheme of approximation can be exploited in many different situations, 
with, of course, differing results depending upon the particular circumstances. For ex- 
ample, in a following work (Eyink (submitted)) we apply our methods to the cascade 
of energy in 2D and obtain results consistent with an inverse cascade. In this case, it 
is a weakly local interaction via the 'skew-strain' which plays the fundamental dynam- 
ical role. We anticipate many other useful applications, such as passive scalar cascades 
and magnetohydrodynamic cascades of energy and magnetic helicity. We hope that our 
method will be useful for all these cases in illuminating the physical mechanisms involved. 
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Appendix A. Convergence of the Gradient Expansion 

We establish here the convergence of the stress approximation r^ n ' m ^ in the limit 
m — > oo, for the space L 1 -norm. The advantage of this norm is that it allows us to 
give the proof by entirely elementary methods (although convergence can doubtless be 
established as well using other L p norms for p > 1). We also give the proof, again for 
simplicity, in infinite volume in d-dimensions without flow boundaries. 

Our argument uses the formula 



_(n,m) 



with 



J d d rG e (r)6u < - n - m ^(r)6u < - n - m \r)-j d d r G e (r)6u^ n ' m \r) J d d r G e (r)Su^ n ' m \r) 

(Al) 



m ^ 



(A 2) 



P =i 



If r has a compactly-supported Fourier transform, then (x) is real-analytic and thus 
5u( n,m \r; x) — » 5u( n \r; x)asnn oo, pointwise in x and also, as we see below, in the 
L 1 norm. Therefore, it is enough to establish absolute integrability and summability: 

r °° i 

d d rG e (r) J2 — ||(r-V) p u (n) (r-V) p 'u (n) ||i < oo (A3) 



p : p'—0 



and 



fd d rG l (r) fd d r'G i (r') -^||(r-V) p u (n) (r'-V/uWyi < oo (A 4) 
J p,p'=0 p ' p ' 

In that case, the integrations and infinite summations commute and 
lim T< n ' m ) = / d d rG t {r)5u^{r)5u^\r) - [ d d r G e (r)6u^ (r) / d d r G e (r)6u^ (r) 

m^oo J J J 

= tW. (A 5) 

Let us establish the bound (A3). By the Cauchy-Schwartz inequality, 

\\(r-V) p u (n ^ (t-V) p V ,i) ||i \r\ p+p ' \\V p u ( ^\\ 2 \\V p ' u^\\ 2 . (A 6) 

Thus, (A3) is implied by 



/ 



d d rG e (r) 



OO 

y | p |p||vp«(")|| 2 



p=0 1 



< OO 



By a similar argument, we see that (A 4) holds, if 



/ 



d d rGi(r) 



OO -. 

y | r |P||vV")|| 2 



Lp=o 



< oo 



(A 7) 



(A 8) 



To proceed we must have an estimate of || V p it(™)|| 2 . This is easy to obtain by going over 
to Fourier transforms using the Planchercl identity: 

\\v p u {n) \\l = J d d k\(ik) p r(e n k)u(k)\ 2 . (A 9) 

Since r has compact support, the integral over k involves only the wavenumbers with 
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\k\ < ci/£ n for some constant c\. Thus, 

||V p M (n) ||2 <; c 2 2 ( Cl /£ n ) 2p J d d k\u{k)\ 2 , (A 10) 

for another constant c 2 = sup fe |-T(fe)|, or 

\\V p u^\\ 2 ^c 2 ( Cl /£ n y\\u\\ 2 . (All) 

With this estimate we get 

oo ^ 

]T-|r|lV^|| 2 ^ C2 exp( Cl |r|/C)||«||2. (A 12) 

p=0 ^' 

Thus, we see that (A 7) and (A 8) hold, if the filter kernel G(r) decays faster than expo- 
nentially in space. For example, G could be Gaussian. Notice that -T(r) can also decay 
very rapidly in space — for example, faster than any inverse power — since the Fourier 
transform r(k) may be both compactly supported and C°°. For any such filter ker- 
nels G and r, we conclude finally that the absolute summability/integrability conditions 
(A 3), (A 4) both hold, and thus the limit relation (A 5) is valid, as claimed. QED. 



Appendix B. Comparison with the Defiltering Approach 

As mentioned in the Introduction, Yeo & Bedford (1988), Leonard (1997), and Carati 
et al. (2001) have developed a somewhat similar gradient expansion for the turbulent 
stress. Here we would like to compare and contrast the two approaches and, in particular, 
indicate our reasons for dissatisfaction with the expansion constructed by those authors. 
Their approach is based upon defiltering, which is the inverse to the filtering operator, 
defined spectrally by 

Q t v(x) = J d d k G{ki)v{k)e lk,x . (B 1) 

Thus, the defiltering operator is given similarly by 

g^v(x) = J d d k [G{k£)]- 1 v{k)e ik ' x . (B2) 

In terms of these operators, Yeo & Bedford (1988), Leonard (1997) and Carati et al. 
(2001) define a tensor-valued functional 

T[v] = Qi{Q^v G^v} -vv. (B 3) 

A little thought shows that if u = Qgu is substituted into this functional, then one 
recovers the turbulent stress as r — T\u]. Hence, this formula provides, seemingly, an 
exact closure of the stress in terms of the filtered velocity u. Furthermore, the functional 
T has a formal gradient expansion: 

TiM =]Tc p , q V^ Vv,, (B4) 

where the summation is over multi-indices p = (p\, ...,p^),q = (qi, qj) with integer 
components. Note that V p = <9f 1 • • • d^ d for the multi-index p. The coefficients in the 
expansion (B4) can be obtained from a generating function (Carati et al. (2001)): 

F[*,1,]= d{ - l{ t + ^ )] =rq,/f. 
GH0)G(-«/>) ^ 
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Here 4> P = 4> V1 ' ' ' 4> P d ■ Thus, the expansion (B4) for r = T\u] seems to yield a closed, 
constitutive formula for the turbulent stress in terms of the gradients of the filtered 
velocity. 

To see the problem with this approach, note that the defiltering operator is not even 
defined if the filter kernel has a compactly supported Fourier transform. In this case, 
all the subscale modes cannot be recovered from knowledge of the filtered velocity u. 
However, this is only an extreme form of a general difficulty. For any filter kernel G, the 
defiltering operator Qj x is unbounded on the natural function spaces for the velocity field, 
such as the LP spaces. This means that defiltering is not defined for every element of those 
spaces, but instead only for a (dense) subspace, called the 'domain' of the operator. The 
natural, maximal domain of the defiltering operator in any of these spaces is the range 
of the corresponding filtering operator on that space, or Dom^ 1 ) = Ran(&). This 
means that, to defiltcr a function in the space, that element must have been obtained by 
filtering another member of the space. However, this is not the case for most functions 
in the space, even infinitely smooth ones. 

For example, consider the most natural space of L 2 or finite-energy fields, and a Gaus- 
sian filter kernel G. Most v e L 2 , even those with very rapidly decaying Fourier coeffi- 
cients v(k), have no Gaussian defilter. Formally, 

[G^vPik) = cxp[+(k£) 2 /2<j 2 ]v(k) (B 5) 

for a Gaussian kernel with Fourier transform G(k) — exp(— fc 2 /2er 2 ). Thus, the Fourier 
coefficients v(k) might decay very rapidly, e.g. exponentially in |fe|, and yet the defiltered 
field has infinite energy or Q7 v ^ L 2 . For such a velocity field, the stress T[v] defined by 
the formula (B3) does not exist. This poses a real difficulty for an LES closure equation 
based upon a constitutive relation t = 7^ m ^[lt] obtained by truncating the expansion 
(B 4) at finite order m. There is nothing to guarantee that the solution u of such a closure 
equation will have Fourier coefficients decaying fast enough for (B 5) to remain in L? . In 
that case, the expansion (B 4) will not converge and there is no reason to expect that the 
solution of the mth-order LES equation will converge in the limit m — > oo to the exact 
filtered velocity, even though the closure then becomes formally 'exact'. Of course, even 
worse, the LES equation may itself be ill-posed and its solution could blow up at finite 
time. The 'exactness' of the closure as m — > oo docs not provide any guarantee of good 
behavior at finite m. 

Such difficulties with defiltering are not unknown in the LES literature. Domaradzki & 
Adams (2002) have reviewed various approaches to subgrid stress modelling by defiltering 
and have pointed out the related fact that defiltering is generally an ill-posed operation 
in function space. In particular, multiplication by the inverse filter transform, as in our 
equation (B5), magnifies the effects of noise and round-off error at high-wavenumbers. 
Thus, the defiltering operation, even when it exists, is not stable to small perturbations 
in the input velocity field. As in any ill-posed problem, various rcgularizations may be 
considered to render it well-posed. Most of the existing approaches have employed an 
approximate, regularized defiltering together with an "eddy-viscosity" or dissipative term 
at high-wavenumbers. See Domaradzki & Adams (2002). Needless to say, our approach 
is quite different. We have constructed an approximation scheme which is proved to 
converge under very modest assumptions but which is not closed. Our present goal is 
to develop a tool to explore the basic physics and not to construct turbulence models 
directly. Our expressions for turbulent stress may be useful in modeling efforts, but we 
must defer to future work the important problem of estimating the unknown subscale 
gradients that appear. 
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Appendix C. CSA Expansion Coefficients for Gaussian Kernel 

We give here the coefficients that appear in the CSA expansion, for the special case of 
an isotropic Gaussian filter: 



Ge(r) 



[2n(<j£) 2 ] d /- 



■ exp 



2{atf 



In general, one would like to have a ~ I, so that this really corresponds to a filter at scale 
i. However, this is often not true; for example, the conventional choice made by Leonard 
(1974) was a 2 = 1/12 (so that the second moment of the Gaussian and box filters would 
agree). In such cases it is better to define the 'shells' in the model formulation as 



S Q = {r : \r\ > 
S k = {r : o-4-i > \r\ > a£ k }, k = 1, 



(CI) 
(C2) 



In this way, increments for separations r w al n are calculated from fields filtered 
at the same length-scale. We must calculate the partial averages of \r\ p , p — 2,4,6... 
with respect to the Gaussian filter, over each of these shells. Introducing a dimcnsionless 
variable p — r/£, these moment averages may be written as cf ] £P for shell S k , k = 
0, 1, ...,n. 

The integrals for coefficients are evaluated by substituting t = p 2 / (2a 2 ) , a 2 dt = 
pdp, and using Sd-i = 2n d l 2 /T(d/2) for the (d— l)-volume of the unit sphere in dimension 
d: 
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for k = 1, ...,n. Here we have introduced the incomplete Gamma functions T(a,x) and 
•y(a,x) = T(a) — F(a,x), as they are defined in the standard literature (Abramowitz & 
Stegun (1964)). 

The asymptotics of these coefficients for large k are obtained using j(a, x) <~ x a /a as 
x -» (e.g. Abramowitz & Stegun (1964), 6.5.4 & 6.5.29). Thus, 

Cf 1 ~ C(d,p,\)\- {d+p)k a p , k->oo, 

with C(d,p,X) = (X d +P - 1)/ [(d + p)T(d/2)2 ( - d - 2 ')/ 2 ] . If we define cjf 1 = \( d +P^C p k \ 
then this new constant becomes independent of k as k 



oo: 



Cf ~C(d,p,\)aP. 



(C5) 



It is not hard to prove that these same asymptotics hold, at least as big-O bounds, for 
much more general filter kernels than Gaussian. 

An interesting application of (C 5) is to establish the order of magnitude of the 'sys- 
tematic' stress term g^'^ that appears in equation (4.13) for T W'( m \ We assume that 
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the velocity field u has the Holder exponent < a < 1. From the definition (4.10) it is 
obvious that Q^'^ m ^ is a sum over p,p' = 0, m of terms proportional to 

cfl p ,£ p k +p '(d p uW)(d p 'uW) (C 6) 

Here d p u^ indicates a pth-order space-derivative of it^l with component indices sup- 
pressed. It is not hard to show (e.g. see Eyink (2005)) that d p v\ k ^ — 0(l'j"~ p ). On the other 

hand, (C5) implies that C p+p , is asymptotically constant as k — > oo. It follows that each 

term p,p' of the Taylor expansion contributes a term for g^'^ that scales as 0(£ 2a ). 
In fact, this is the correct order of magnitude for the contribution to the stress from 
length-scale £k (Eyink (2005)). It is worth pointing out that the term u}^' 1 ^ 1 u} k ''^ m ' 1 
in (4.13) would likewise scale as 0(£^£^,), if the correction factor Nk had been used in the 

definition (4.12) rather than N]J 2 . Indeed, so defined, u'} k ^ m " l u'} k ''^ m ' 1 is from (4.10) 

a sum over p,p' of terms proportional to cjf'Cy ^ £ p .£ p ,(d p u^)(d p ' u^ k '^). In that case, 
the previous argument carries through. It might seem that this alternative definition of 
u '[k],{ m ) th ere f ore S ome merit, except that it ignores the cancellations that we expect 
to occur in the integrals over volume. 



Appendix D. Generalized Betchov Relation in 3D 

We here give briefly the proof of the generalized Betchov relation (5.12). Let us suppose 
that a, b, c are three solenoidal (divergence-free) vector fields in 3D. Then, it is trivial to 
verify from the product rule of differentiation that 

'i... 1 ',. >,<'>,., + CkjbjjOi^ = d l (a l jbj, k Ck) ~ dkiaijbj^Ck) + dj(ai^bj^Ck) (D 1) 

where we use the notation aij = dai/dxj, etc. Because the righthand side is a total 
space-derivative, the ensemble-average of the lefthand side is zero if one assumes space- 
homogeneity: 

(dijbj^Ckj + Ck,jbj t i<ii t k} = 0. (D 2) 

This same relation holds for space-averages, if boundary conditions on a, 6, c are such 
that boundary terms from integration by parts can be ignored. With either of these 
assumptions, let us then apply (D 2) for b = u and a = c = u^ n \ This gives 

(Tr[D{D {n) ) 2 }) =0, (D3) 

with Dij = Uij the deformation matrix associated to a velocity field u. 
Now, using D-™^ = — (l/2)ejjfco>(™) [equation (5.3)] gives 

(D {n) f = (S (n) ) 2 - i(/|u/ n) | 2 - u) {n) u) {n) ) 

+ i(u/")xS (rl) +S (n) Xt>)). (D4) 

Note that the matrices on the first line of the righthand side are symmetric and the 
matrix on the second line is antisymmetric. (/ is the identity matrix.) Thus, if one uses 
the relation analogous to (5.3) for D^, then the only contribution to the trace of D 
with (D^) 2 is from the trace of Sij with the first line in (D4) and from the trace 
of (l/2)eijfc5Jfc with the second line. Using tracelessness of the strain matrix S and the 
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identity e^eum = Sji5 km - S jm S k i gives 

Tr[D(D (n) ) 2 ] = Tr (s(S (n) ) 2 ) + J(u;< n >) T 3(a;< n >) + ^(ZJ) T S (n) (u: (n) ) . (D 5) 

Substituting (D5) into (D 3) gives (5.12). QED. 
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